Non-Arrhenius Behavior of Surface Diffusion Near a Phase Transition Boundary 
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We study the non-Arrhenius behavior of surface diffusion near the second-order phase transition 
boundary of an adsorbate layer. In contrast to expectations based on macroscopic thermodynamic 
effects, we show that this behavior can be related to the average microscopic jump rate which in 
turn is determined by the waiting-time distribution W(t) of single-particle jumps at short times. 
At long times, W(t) yields a barrier that corresponds to the rate-limiting step in diffusion. The 
microscopic information in W(t) should be accessible by STM measurements. 



PACS numbers: 68.35.Fx, 82.20.Pm 
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The migration of atoms and molecules is one of the 
most important processes taking place on solid surfaces. 
It appears in many phenomena such as catalytic reac- 
tions and surface growth that are important for practical 
applications p]. In most experimental and theoretical 
studies of the surface diffusion constant D, its temper- 
ature dependence is analyzed through an assumed Ar- 
rhcnius form, where D is written as a product of an 
entropic prefactor Dq and a term exp(— /UbT) de- 
scribing thermally activated jumps over an energy barrier 
E®. Although the Arrhenius form can be derived from 
microscopic considerations in some special cases HH, a 
rigorous justification for its use in interacting systems 
at finite coverages is not available. Further, even in the 
cases where D appears to have an Arrhenius tempera- 
ture dependence over a finite temperature range, its mi- 
croscopic interpretation may not always be clear. This 
is because for an interacting system, there may be many 
microscopic activation barriers. Thus the value of the 
measured effective diffusion barrier E^ must result from 
some complex average of all of them, and does not refer 
to any microscopic process in particular H . 

In fact, the values for Dq and can be strongly 
temperature-dependent indicating a region of non- 
Arrhenius behavior. This becomes especially pronounced 
near surface phase transition boundaries, where rapid 
variations of D have been observed in experiments 
and computer simulations Such rapid changes are 

often accompanied by the well-known "compensation" ef- 
fect Q , where an apparent increase in E® is compensated 
by an increase in the prefactor Dq H- However, in most 
cases the underlying reasons for non-Arrhenius behavior 
are not understood. It is the purpose of the present work 
to study these issues near a second-order phase transi- 
tion in a surface adsorbate layer. We show that in con- 
trast to the common folklore that an anomalous temper- 
ature dependence in D near T c would be predominantly 
due to non-local thermodynamic effects, it can be ex- 



plained by the microscopic single-particle jump rate T. 
This quantity is determined by the short-time behavior 
of the waiting-time distribution W(t) for single-particle 
jumps. Moreover, we show that for long times, W(t) 
yields an effective activation barrier that corresponds to 
the rate-limiting step in diffusion. Thus W(t) provides 
a connection between microscopic and macroscopic as- 
pects of diffusion. Further, it is experimentally available 
through e.g. STM measurements ||. 

In this Letter, we have carried out Monte Carlo (MC) 
simulations for a model of oxygen on the W(110) surface 
|]l0l , [ri"l| . In this system, the substrate remains unrecon- 
structed the oxygen atoms have well-defined adsorp- 
tion sites [|13fl , and desorption of oxygen occurs only at 
temperatures 1600 K or above [Q. Therefore, this sys- 
tem is very suitable for simulation studies using a lattice- 
gas description. We use the lattice-gas model constructed 
by Sahu et al. |ll|] to describe the main features of the 
phase diagram. The Hamiltonian includes pair interac- 
tions up to fifth nearest neighbors and some three-body 
interactions pi], the attractive ones being the dominant. 
We concentrate on results for the coverage 8 — 0.45 over a 
wide temperature range. For this coverage at a low tem- 
perature, the adlayer is in the ordered p(2 x 1) phase, 
while at T c « 710 K it undergoes a second-order tran- 
sition jll|] to a disordered phase |H. For details of the 
model and MC simulations, see Refs. [ ^0p^ ]. 

Our simulation results for the tracer and collective dif- 
fusion coefficients Dt and Dq (for definitions see e.g. 
Ref. pj), respectively, are given in Fig. |l]. We first note 
that their qualitative behavior is similar and that the 
effective diffusion barrier E® defined as 

D _ d(\ogD) 



d(l/k B T) 



(1) 



is approximately constant at low and high temperatures 
away from T c . This implies that the diffusion constants 
obey simple Arrhenius behavior. Near T c , however, the 
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temperature dependence of the diffusion constants is 
strongly non-Arrhenius. 

In the Green-Kubo formalism the expression for 
Dc contains a thermodynamic factor £ inversely propor- 
tional to the compressibility, which is governed by the 
global number fluctuations of the adlayer. It is often as- 
sumed that a possible non-Arrhenius behavior of Dc near 
the phase transition boundary is predominantly due to 
the critical behaviour of this factor . We show below 
that this is not the case here: The non-Arrhenius behav- 
ior of both Dt and Dc has a dynamic origin and can 
be traced back to the temperature dependence of the 
local single-particle jump rate T. To demonstrate this, 
we show in Fig. [I] the temperature dependence of £ and 
r as well as that of Dt and Dc- It can be seen from 
Fig. [j] that while the diffusion constants have a turning 
point and sharp temperature variations close to T Cl the 
thermodynamic factor £ has only a relatively weak tem- 
perature dependence in this region and cannot account 
for the non-Arrhenius behavior of the diffusion constants. 
On the other hand, the single-particle jump rate T has 
exactly the same behavior near T c as Dt and Dc ■ These 
observations can be understood theoretically within the 
dynamical mean field theory (DMF) @, which yields 
Dt oc r and Dc oc £T. We can conclude that the strong 
temperature dependence of both Dt and Dc near T c is 
indeed of the same dynamic origin, coming from the av- 
erage single-particle transition rate T. 

We next focus on the effective diffusion barrier as 
extracted from Eq. ([!]) for Dt- As shown by squares in 
Fig. ||, has a sharp peak near T c . This peak in E® 
is accompanied by a strong increase in the value of the 
corresponding prefactor Dq shown in the inset of Fig. |^. 
This is yet another example of the well-known compen- 
sation effect . Here the compensation simply results 
from the fact that when the temperature dependence is 
non-Arrhenius, there is no unique way of separating the 
prefactor and the barrier contributions. Since the tem- 
perature dependence of the diffusion constant itself near 
T c is smooth and nonsingular, any dramatic change in 
the temperature dependence of the effective barrier E® 
must be followed by a corresponding change in the effec- 
tive prefactor Dq. We note that the same phenomenon 
occurs for collective diffusion as well. 

To understand the observed strong temperature vari- 
ation of E® near T c , we need to consider the energetics 
of the microscopic jump processes which determine the 
average jump rate T. At finite coverages, there is a very 
complex distribution P(E a ) for the instantaneous activa- 
tion barriers E a Jl7[| which an adatom needs to overcome 
in a jump attempt from one configuration to another. At 
high T, P(E a ) is strongly peaked at small values of E a , 
while at low temperatures the situation is completely the 
opposite |l^] . The change in the distribution takes place 
around T c , thus characterizing the ordering of the adlayer 
as the temperature is decreased below T c . This change in 



turn results in a strong temperature dependence of the 
average transition rate T around T c , as shown in Fig. |]. 
We point out that the instantaneous activation barriers 
E a cannot explain the peak of the effective barrier E® in 
Fig. |^, since the largest value of E a in our model system 
is only about 0.4 eV |T^]. Thus, the peak does not re- 
fer to any microscopic rate-limiting process. Instead, it 
arises from an entropic contribution |L9| to T which has 
a strong temperature dependence in the vicinity of T c . 

To gain more insight into the microscopic dynamical 
processes and the anomalous temperature dependence 
near T c , we next introduce the waiting-time distribution 
W(t) of single-particle jumps ^o). Suppose a single par- 
ticle (in the presence of other particles) had performed 
its last transition at time t — 0. Then W(t) is the prob- 
ability density that the particle in question performs its 
next transition at time t after it remained still until t. 
Here the most practical definition of "time" in the MC 
simulations is to consider the time scale as the number 
of jump attempts of the particle, denoted by n. Then 
the waiting-time distribution is simply W(n). This pro- 
vides a direct connection with the dynamic jump rate T 
discussed above via 

^ oo 
n— 1 

where (n) is the average waiting-time of the particle. 

At very long times, we expect W(t) to decay as W(t) ~ 
exp(— i/r). Here the characteristic time r describes the 
longest time scale among the various microscopic pro- 
cesses, which constitutes the rate-limiting factor for mass 
transport. This expected exponential decay at long times 
is indeed observed for our model system, as demonstrated 
in Fig. |j. We can then define an effective activation 
barrier E^ via r by considering the jump probability 
p = 1/t = po cxp(— EY /ksT). As shown by circles 
in Fig. ||, the activation barrier E^ extracted from the 
asymptotic region of W(t) decreases monotonically with 
increasing temperature, and agrees with the effective dif- 
fusion barrier E® extracted from an Arrhcnius analysis 
of Dt far from T c . Additional studies in our model sys- 
tem indicate that the value of E^ is closely related 
to the instantaneous activation barrier characterizing the 
dominant microscopic processes. In our model the micro- 
scopic barriers have a maximum value of about 0.4 eV 
and thus the barrier Ejf does not have the sharp peak 
displayed by the effective diffusion barrier E%. 

It turns out that the temperature dependence of the 
barrier E® results mainly from the short-time behavior 
of W(n). This is demonstrated by dividing the sum in 
Eq. (g) into two parts, the first of which is the short- 
time contribution (n)s — Yln=i n W(n). This quantity 
accounts for the contribution up to a crossover time n co , 
which separates the short-time regime from the asymp- 
totic exponential decay. What remains is the long-time 
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contribution (n)i, = (n) — (n)s- As expected, from Fig. |] 
we observe that the short-time regime gives the domi- 
nant contribution to T. Further, the short-time regime 
of W(n) is strongly affected by the critical fluctuations, 
being mainly responsible for the anomalous temperature 
dependence of the diffusion constants near T c . 

To summarize, within the present model of 0/W(f f 0), 
the non-Arrhenius behavior near T c was found to have 
a mainly dynamic origin, reflecting the dependence of 
the single-particle jump rate T on the critical fluctua- 
tions close to T c . Surprisingly, in our studies the ther- 
modynamic factor gives only a minor contribution to the 
temperature dependence of the collective diffusion con- 
stant Dq, and the anomalous temperature dependence 
for both Dt and Dc results from the dynamic factor 
r. We find that the single-particle waiting-time distribu- 
tion W(t) gives the most detailed picture of the micro- 
scopic processes. It has been recently demonstrated by 
Swartzentruber that this distribution function in the 
presence of several different microscopic activation bar- 
riers can indeed be measured using the STM. From the 
long-time tail of W(t) one can obtain information on the 
energetics of the rate-limiting processes of diffusion in the 
form of an effective activation barrier E^ . On the other 
hand, the temperature variation of the effective diffusion 
barrier directly reflects that of the microscopic jump 
rate T, and depends not only on the long-time tail of 
W(t) but also on its short-time behavior. In the short- 
time regime near T c , W(t) is strongly affected by the crit- 
ical fluctuations. The fluctuations there lead to a strong 
temperature dependence of the transition entropy and an 
additional contribution to the effective barrier E^. 
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FIG. 1. Results for Dt and Dc as a typical Arrhenius plot 
in the O/W(110) system at 6 = 0.45. We also show the behav- 
ior of the average transition rate V and the thermodynamic 
factor £. The quantities Dt, Dc, and T have been made di- 
mensionless by expressing them in units of a 2 /to, a 2 /to, and 
I /to, respectively, where a is the lattice constant and to is 
one Monte Carlo step per particle. The value of £, which is a 
dimensionless quantity, has been scaled by a factor of seven 
to make the graphic representation more readable while other 
quantities are expressed directly in the units given above. The 
critical temperature of the order-disorder phase transition is 
denoted by T c and a dotted line. 
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FIG. 2. Results for the effective activation barriers. The 
squares denote results based on the Arrhenius form (see 
Eq. (^)) for tracer diffusion, while open circles represent the 
data based on the tail of W(n). Behavior of the prefactor Do 
is illustrated in the inset. The critical temperature is denoted 
by a dotted line. 
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FIG. 3. An example of a waiting-time distribution W(n) 
at a temperature of 0.774 T c showing an exponential decay at 
long times. For clarity, only some of the data points are shown 
here. The full curve is an exponential fit to the tail of W(n). 
The approximate crossover time n co for the crossover from the 
small-time regime to the asymptotic long-time regime is 
indicated by an arrow. 
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FIG. 4. Comparison of the short-time contribution (n)s 
and the long-time contribution (n)L to the average wait- 
ing-time (n). The slight increase of (ti}l at small T is due 
to n co whose value is difficult to determine accurately at very 
low temperatures. The quantities (n)s, (n)h, and (n) are all 
expressed in units of one Monte Carlo step per particle. 
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